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METHOD FOR RECONSTRUCTING 
THE TOPOLOGY OF A POLYGONAL SOUP 

. _ Technical Field 

The present invention relates to 
5 visualization tools for interactive simulations with 
complex vehicle models. 

Background Of The Invention 

To meet the need for interactive simulations 
with complex vehicle models, a generation of 

10 visualization tools have been developed. However, most 
of these tools have limitations inherent in their data 
structure. In this regard, most visualization tools 
have adopted a polygonal soup, that is, an unstructured 
collection of polygons, as their standard data 

15 representation. However, a variety of existing 

applications, such as CFD (Computational Fluid 
Dynamics) , FEA (Finite Element Analysis) , and assembly 
simulation, require more than a cloud of points or 
triangles. 

2 0 Thus, a need exists for a method of 

reconstructing the topographical information 
(relationships between polygons) of a polygonal soup in 
order to overcome the known problems with visualization 
tools. 

2 5 Summary of the Invention 

The object of the present invention is to 
provide an improved visualization tool for interactive 



simulation with complex vehicle models. It is another 
object of the present invention to provide a 
visualization tool that comprises more than a cloud of 
points or triangles. 
5 . . - in accordance with the present invention, an 

algorithm is provided which automatically reconstructs 
topological information for a given mesh and then 
alters the mesh by introducing, deleting, or splitting 
existing polygons when needed. 

10 In accordance with the present invention, a 

new OctTree space decomposition is used to achieve a 
log2 -complexity search which locates the closest vertex 
in the polygonal soup to a given point in space . A 
technique with linear complexity is then used to locate 

15 all of the triangles connected to that vertex. A 
technique with linear complexity is also used to find 
all triangles connected to a given triangle. The 
triangles are then split to enforce conductivity. 

These and other objects, features, and 

2 0 advantages of the present invention will become 
apparent from the following detailed description of the 
invention when viewed in accordance with the 
accompanying drawings and appended claims. 

Brief Description Of The Drawings 

2 5 FIGURE 1 depicts an exemplary cube 

representation for modeling tessellated meshes; 

FIGURE 2 sets forth index representations of 
the triangles and vertices of the cube depicted in 
Figure 1; 
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FIGURE 3 depicts a dynamic vector vertex- 
neighbor table; 

FIGURE 4 depicts tables generated from the 
vertex-neighbor table; 
5 - • ~ FIGURE 5 depicts a dynamic vector edge- 

neighbor table; 

FIGURE 6 illustrates an object compound of 
unconnected strips of connected triangles ; 

FIGURE 7 depicts a generalized example of a 
10 step in synchronizing triangle strips; 

FIGURE 8 depicts the formation of new 

triangles; 

FIGURE 9 illustrates the exemplary cube as 

modified; 

15 FIGURE 10 illustrates another step in the 

sorting process; 

FIGURE 11 illustrates another step in the 
sorting process; 

FIGURE 12 depicts another sorting example; 
20 FIGURE 13 depicts still another sorting 

example; 

FIGURE 14 illustrates an OctTree example; and 
FIGURE 15 depicts a further sorting example. 

Description Of The Preferred Embodiment(s) 

2 5 As indicated, there is a need for interactive 

simulation with complex systems, such as airplanes and 
cars, and in fact visualization tools have been 
developed which are capable of rendering large models 
in real-time. However, in developing the 3-D rendering 
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tools, insufficient attention has been paid to 
architecting systems capable of evolving into effective 
modeling tools. One of the critical limitations is 
inherent in the data structure. Almost all 

5 visualization tools have adopted polygonal soups as 
their standard data representation. However, a variety 
of existing applications, ranging from CFD, FEA and 
assembling simulations, require more than a cloud of 
points or triangles. 

10 Currently, different data representations are 

used to model tessellated meshes. Figure 1 illustrates 
one example represented by a cube 10. The cube is 
modeled as a set of eight triangles, two per side, and 
eight vertices, V 0 through V 7 . 

15 In the STL (stereo lithography) format, the 

model is described in either ASCII or binary format, as 
a list of triangular elements. The actual coordinates 
of the three vertices of each of the triangles are 
expressly stored and are shown in the lefthand table 

2 0 set forth in Figure 2. The actual coordinates of the 

three vertices are provided on the left side of the 
table in Figure 2, and the cube is represented as a 
list of 3 X 12 vertices. In this regard, other 
representations, such as Virtual Reality Modeling 
25 Language (VRML) , OpenGL, Optimizer, DirectX, and 
DirectModel, offer other efficient, indexed 
representations . 

Figure 2 shows two separate lists,- a vertex 
list with the coordinates of all vertices, and a 

3 0 triangle list where the indexes of the vertices for. 

each triangle, relative to the vertex list, are stored. 
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The STL representation is, in fact, a true 
polygonal soup, that is, the geometry is a set of 
unconnected triangles. The index representation is 
better. In this regard, the same vertex is shared 
5 among ►multiple triangles and provides some topological 
information. However, the facets are usually connected 
in strips or fans for graphics performance, without a 
collective topological map and . with multiple 
representations of the same vertex on different strips. 

10 This causes small gaps in the geometry to appear 
because of numerical approximations. 

Real-time tools, such as collision detection 
tracking algorithms, take advantage of spatial 
coherency between successive time samples. This 

15 requires information about vertex and edge 
connectivity, such as the neighbors of a given triangle 
and of a given vertex. For example, in the case of the 
cube shown in Figure 1, the triangle T 4 is connected to 
triangles T x , T 2 and T 5 , and the vertex V 2 is shared 

2 0 among the triangles T 1# T 2 , and T 4 . 

The present inventive method basically 
consists, of three steps or phases. In the first phase, 
the vertex and edge conductivity information is 
developed starting from available information. The 

2 5 second and third phases deal with imperfect meshes. 

During the second phase, vertex duplicates are 
discovered and eliminated. During the third phase, the 
model is re-meshed to realign strips of triangles that 
do not share common vertices. 

3 0 The initial phase is divided into three 

additional steps. In the first step, the index 
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representation of Figure 2 is generated. An initial 
attempt of reducing duplicated vertices is performed at 
this stage. In this regard, the first step operation 
is only required when the geometry is described as an 
5 explicit set of triangles, such as STL. In the second 
and third steps, the vertex-neighbors, and edge- 
neighbor tables are developed. 

The vertex-neighbors table is a dynamic 
vector where, for each vertex, the list of all the 

10 connected triangles is immediately available, as is 
shown in Figure 3 . The algorithm sequentially browses 
the triangle list. For each triangle, the table listed 
in Figure 2 contains the indexes of the three vertices. 
An identifier, pointing to the triangle sequential 

15 number in the triangle-list is added to the index- 
neighbor list, in correspondence to each vertex. For 
example, in the case of the cube shown in Figure 1, 
this method would start by looking at the triangle T 0 
with vertices V 0 , V x , and V 3 . Table A in Figure 4 is 

2 0 then generated. 

Proceeding with T x , the additional entries are 
added to Table B. Thereafter, the complete vertex- 
neighbor Table C is developed. Since the triangle- 
index table is traversed only one time, this creates an 
25 efficient, linear-complexity algorithm. In this 

regard, the number of computations is n t -x3, where n t is 
the number of triangles. 

The edge ^neighbors table is a vector of lists 
with an entry for each triangle. Each list has three 

3 0 elements, one for each edge of the given triangle, "V X1 - 

v i2' v i2" v i3/ and v i3" v io/ respectively, depending on the 
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position of the element in the list. These elements 
can be equal to either -1, to represent an unconnected 
edge, or to the sequential number of the connected 
triangle to the selected edge, as can be seen in Figure 
5, 5. -Thereafter, the algorithm sequentially browses the 
triangle list, starting with the first element, and 
determines if there are other triangles that share any 
of its three edges. For example, with the cube 
illustrated in Figure 1, the first edge of the 

10 first triangle T 0 can be viewed as follows: 



Triangle 




v i2 


y i3 


T ° 


v 0 




v 3 



(1) 



Two triangles are connected at an edge if 
they share the same two vertices making up the edge 
15 extremes, in this case V 0 and V x . Then, a list of 
potential candidates is found looking for all triangles 
connected to V 0 that are not T 0 in the vert ex -neighbors 
table. This is shown as follows: 



Vertex 




T i2 








v 0 


T 0 


T 8 


T 9 


T 10 





(2) 



20 

Since there are now four triangles, Tg-T^, it 
is necessary to look at the triangle table to see if 
any of the candidates contains the vertex V x as well: 
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Tri angl e 






v i3 


T 8 




v 0 


v 5 


T 9 


v 4 


v s 


v 0 


Tio 


v 3 


v 7 


V 0 




v 4 




v 7 



(3) 



In this particular example, the first edges 
of T 0 and the first edge of T 8 share the same vertices 
V x and V 0 . The two triangles are connected, and the 
5 edge -neighbor table is then updated as shown below: 



Triangle 


T ±1 


T i2 


T i3 






T e 
















(4) 


T 8 


T 0 








As a result 


, after 


browsing 


all of the 


edges 



of all of the triangles, the following table is 
10 completed: 
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Tri angl e 




T 12 




T 0 


T s 


Tx 


T 10 


Tx 


T 2 


T 0 


T 4 


~ • T 2 


T 4 


T 3 


Tx 


T 3 


T 10 


T 2 


T 7 


T < 


T x 


T 5 


T 2 


T 5 


T 7 


T 4 


T a 


T 6 


Txx 


T 7 


T 9 


T 7 


T s 


T 6 


T 3 


T B 


T 0 


T 9 


T 5 


T 9 




T 8 


T u 


T 10 


T 3 


Txx 


T 0 




T 9 


T 10 


T 6 


As shown in 


Table (5) , 


it is 


noted that 



are not any unconnected edges and, therefore, this 
5 particular geometry is a manifold body. 

The total number of computations depends on 
the average number of triangles connected to a given 
vertex. The triangle/vertex ratio n t/v ranges from a 
value of three in the case of a tetrahedron, to 4 . 5 for 
10 the cube described above, to about 5.0 for typical 
automotive models. This means that the number of 
computations to find the triangle connected to a given 
edge is: 

. (n t/v -l) -3 (6) 
15 Since each triangle has three edges, and there are n t 
triangles, then the total number of computations is: 

n t - <n t/v -l) -9 « n t -36 (7) 
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This provides an efficient, linear-complexity 
algorithm. 

The algorithms used to reconstruct the 
topographical information of a tessellated mesh, as 
5 discussed above, rely on two assumptions. The first 
assumption is that duplicate vertices have the same 
coordinates. This requirement is a foundation of a 
variety of search engines, usually based on balanced 
trees or skip-lists, used to detect duplicate vertices. 

10 The second assumption is that the tessellation is built 
without gaps. This translates in having strips or fans 
of triangles all topologically connected at the same 
set of vertices. However, in particular with VRML and 
OpenGL files, the objects are composed of unconnected 

15 strips of connected triangles, as shown in Figure 6. 
This prevents building of edge connectivity since two 
triangles are considered connected at an edge only if 
they share the entire edge. 

As a result, existing visualization tools for 

20 complex systems have a plurality of gaps, poor 
alignment, and degenerate meshes. Also, duplicate 
vertices and the synchronization of triangle strips are 
not sufficiently taken into account. 

In the present invention, these problems with 

25 existing systems are considered and taken into account. 
Duplicate vertices are eliminated and triangle strips 
are synchronized. The algorithm for accomplishing this 
is described below. 

If the vertex V 0 is duplicated and the copies 

3 0 replaced by the vertex V ob in the cube model described 
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above to describe the triangles T 10 and T L1 , the new 
vertex and triangle list are then modified as follows: 



Vertex 


X 


Y 


Z 


Triangle 




V i2 V i3 


V 0 - - 


• -1 


-1 


+ 1 








V ob 


-1 


-1 


+ 1 


Tio 


v 3 


V 7 V ob 










Tn 


v 4 


V ob V 7 



(8) 



10 Due to the new duplicate vertex, the vertex-neighbors 
and edge-neighbors tables are also changed to the 
following: 
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Triangle 


Tii 


T i2 


Tit 


Vertex 


Tn 


Ti2 -~ — »• 


To 


T 8 


Ti 


-1 




V 0 


To 


T 8 T 9 














V ob 


■Tio 


Tn 


T 9 


T s 


T 9 


-1 












Tio 


T 3 


Tn 


-1 












Tn 


-1 


Tio 


T 6 



(9) 



20 

As a result, all four triangles, T 0 , T 9 , T 10p and T 1X , now 
each have an unconnected edge. 

In order to remove duplicate vertices, the 
algorithm first browses the edge table searching for 

25 unconnected edges. Once an edge is found, duplicate 
copies of either one of the two vertices at the ends of 
the segment are looked for, by searching in the vertex 
list for the closest vertex. In Table 8 above, 
vertices V 0 and V ob have the same coordinates, and 

3 0 therefore have minimum distance, equal to zero. 
However, in other cases, the distance is small, but 
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significant. It is necessary to make sure that the mesh 
will not degenerate once a vertex is replaced with 
another one . 

Once a duplicate vertex has been found, the 
5 procedure is as follows: (a) in the triangle index 
table, all references to the duplicate vertices are 
replaced with a reference to the original vertex; (b) 
the entry in the vertex table for the duplicate vertex 
is removed and the appropriate connected triangles are 

10 added to the entry list of the original vertex; and; 
and (c) the edge table for all triangles connected to 
"original" vertices are rebuilt. When applied to the 
copy of duplicated vertices, V 0 and V ob , this sequence 
restores the original vertex, triangle, vertex- 

15 neighbors and edge -neighbors tables. 

The problem with synchronizing triangle 
strips, such as those depicted in Figure 6, can be 
generalized as shown in Figure 7 as a vertex falling 
inside another triangle edge. The introduction of two 

2 0 new triangles, namely T la and T lb as well as a new vertex 
V 8 , changes the vertex and triangle table as follows: 





vertex 


X 


Y Z 


25 


v 9 


0 


0 +1 



Triangle 




Vi2 


Vu 


T la 


v 2 


Vg 


Vi 


Tib 


v 2 


v 3 


Va 



(10) 



30 



As evident from the new edge-neighbor table, there are 
three edges that are not connected: 
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Triangle 


Tu 


T i2 


T i3 




T 8 


-1 


T 10 




T lb 


-1 


T 4 


T lb 


T 2 


-1 





There are different alternatives for building 
a topological connection. First, an artificial vertex 
could be introduced in the triangle T 0 . However, this 
5 would require dealing with polygons having different 
numbers of sides since with a new vertex, T 0/ would be 
a quadrilateral. Secondly, the vertex V a could be moved 
to overlap either one of the two vertices V x or V 3 , and 
then be eliminated. However, this approach would only 

10 modify the geometry, with the potential loss of 
important details, if the triangles connected to V 8 do 
not share the same normal. Thirdly, the triangle T 0 
could be split into two new triangles, replacing the 
edge V X V 3 with V^g and V 3 V 8 , respectively. 

15 The present invention implements, the latter 

alternative. As in the method dealing with duplicated 
edges, the algorithm browses the edge table looking for 
unconnected edges. Once an edge is found, it then 
looks for any vertex that falls very close to the edge 

2 0 itself, but not on one of its vertices. In the example 
shown in Figure 7, the first unconnected edge that is 



Triangle 


Tn Ti 2 Ti 3 


Tri angle 


Vii v i2 v i3 

To 


T 8 -1 Tio 
Tib -1 T 4 
T 2 -1 T 10 


To 


V 0 Vx_ V 3 

•L la 
T lb 


25 



(12) 
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encountered is the second edge V 1 V 3 of the triangle T 0 : 
Searching in the vertex list, the closest vertex to 
this segment is V 8 . Since the distance to the segment 
is equal " to zero, and the distance to both segment 

5 extremes - is equal to 4l , the vertex is recognized to 
fall on an inside edge. Once the closest vertex is 
found, it is necessary to search in the vertex- 
neighbors table to see if there is a triangle connected 
to it which has an unconnected edge parallel to the 
10 original edge. Moreover, the triangle must not overlap 
the original triangle. In the "cube" example discussed 
above, there are two triangles connected to V 8 : 



Vertex 




'J. ' ^ g ••• *•• ••• 


v 8 


T la . 





(13) 



15 Once a triangle, edge, and vertex have been 

identified, the geometry is then re-meshed by splitting 
the triangle. First, the triangle table is updated to 
reflect the fact that the original triangle is split 
into two new elements. In general, the original 

2 0 triangle is kept with one of its vertex indexes 
replaced and a new triangle is added. The triangle 
table will then change as follows: 



25 



Triangle 


1'xj.ajjyle 

"ii V i2 V i3 


Vil v i2 V i3 


T'o 


Vo V L V 3 io 


V 0 Vx V 8 
Vo V 8 V 3 


i l2 



(14) 
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In the generic case, the algorithm must take 
into account which edge the unconnected vertex falls 
on, as shown in Figure 8. The generic rules to update 
the triangle table are then summarized as follows: 
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Before 


After 




Triangle 


Vii 


V i2 


v i3 


v u 


v i2 


V i3 


edge 


Torig 


Vii 


V i2 


V i3 


Vii 


Vnew 


v i3 


Vu -V ia 


T n ew 








V n ew 


v i2 


V i3 


edge 


Torig 


Vii 


v i2 


V i3 


Vii 


v i2 


Vnew 


V i2 -Vi3 


Tnew 








Vii 


V n ew 


V i3 


edge 


Torig 


Vu 


v i2 


V i3 


V netf 


v i2 


v i3 


Via-Vu 


Tnew 








Vii 


V i2 


Vnew 



(15) 



Thereafter, the vertex- neighbor table, is 
15 updated. In the case of the cube shown in Figure 7, 
this would change as follows: 



20 



25 



30 



Vertex 




T i2 










V 0 


To 


T 8 


T 9 


Tio 


Tn 




Vi 


To 


Tia 


T 4 


T s 


T 8 




v 3 


To 


Tib 


T 2 " 


T 3 


' Tio 




v 8 


Tia 


Tib 












Vertex 


Tn 


T i2 










V 0 


To 


T 8 


T 9 


Tio 


Tu 


T12 


Vi 


To 


Tia 


T 4 


T 5 


T 8 




v 3 


T lb 


T 2 


T 3 


T 10 


T12 




V 9 


T la 


Tib 


T o 


T12 







(16) 
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In the general case, the following table 
summarizes the rules needed for determining the new 
table, depending on which of the three angles is 
affected by the new vertex: 





Be fore 


After 




Vertex 


Tii 


Tii 


edge 


Vii 








Vil-Via 


v i2 


... Torig 


T ne w 






V i3 
Vnew 




Tnew 
Torig 


Tnew 


edge 


Vii 




Tnew 




v i2 -v i3 


v i2 










V i3 
V ne w 


Torig 


Tnew 
Torig 


Tnew 


edge 


Vii 


T 0 rig 


Tnew 




V 13 -Vii 


V i2 




- . Tnew 






V i3 
V n ew 




Torig 


Tnew 



(17) 
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The last step is to update the edge -neighbors 
table in correspondence of the three vertices of the 
original triangle V 11# V i2 and V i3 as well as the vertex 
projecting in the middle of the edge V new , using the 
algorithm described above in equations (6) and (7) . In 
the case of the present example, the edge -neighbors 
table is updated as follows: 



Triangle 




T i2 


T i3 


T 0 


■T 8 


T la 


T 12 


T la 


Ti b 


T 0 


T 4 


T lb 


T 2 


T 12 


.T la 



(18) 



16 



10 



The geometry is now manifold as shown in Figure 9. 

An efficient closest-vertex search is 
utilized in accordance with the present invention. If 
the entire geometry is searched, the number of 
calculations required to eliminate duplicate vertices 
and split triangles in order to obtain a connected, 
closed manifold is a function of the square of the 
number of vertices and shown by the following equation: 
a a n. 2 a n e/v 2 _ a 3 n t/v 2 

where n e is the number of edges, n„ is the number of 

vertices, is the edge/vertex ratio, ^y^yy is the 

triangle/vertex ratio, and a is the percentage of edges 
that are not connected and require some action. Since 

15 a typical tessellated model of a single vehicle 
component usually consists of hundreds of thousands of 
triangles, the number of computations quickly grows and 
becomes unmanageable when dealing with complex 
assemblies. Different solutions are available. Most 

2 0 of these are based on the same philosophy where all of 
the elements of the model are assigned to some bounding 
volumes of simple geometry. In order to reduce the 
number of. computations from a linear to a logarithmic 
dimension, the hierarchy of bounding volumes is often 

2 5 used. Some implementations subdivide the space in 
spherical -containers, while others make use of 
variable-side bounding boxes. 

With the present invention, the space is 
decomposed using an OctTree structure. This is a data- 
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structure, similar to a voxel map, and is in current 
use in computer-graphics. An example illustrates its 
concept. If the parent mode of the OctTree is a box 
enclosing the entire workspace, vertices are assigned 
5 to -it.- - As soon as there are more than a constant, 
assigned number of vertices in the parent node, for 
example 4, the box is split into a set of 8 identical 
smaller boxes, where the intersection is the empty set 
and the union is the original box. The elements 

10 previously assigned to the original box are then 
reassigned to the smaller containers. This is shown in 
Figure 10. As the sorting continues and new vertices 
are assigned to the proper boxes, the sub-box that is 
filled is once more subdivided into eight smaller 

15 volumes and so on. This is shown in Figure 11. 

The data structure is bi-directional. Each 
node has a pointer linked to either a list of vertices 
or. a list of sub-boxes. On the other hand, each vertex 
has embedded in its coordinate the address of the box 

20 containing it. Since the boxes are always split and 
each side is dissected into two equal segments, it is 
sufficient to find the address of the box to normalize 
the coordinates by the size of the workspace and then 
multiple the result by 2 n , where n is the maximum number 

2 5 of layers, and round off to the nearest integer. For 
example, given the point shown in Figure 12, the 
address in an n=4 layers OctTree is as follows: 



int 



)4 



V "^rnax *min / 



,int 



y - " y ™ . -2 4 => 001 1,0001 (20) 

V -^max ^min / 



These indexes can be used to find the 
3 0 containing box. If the parent node points to a list of 
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vertices, then the search is completed. Otherwise, the 
Most Significant Bits (MSB) defines the address of the 
sub-box which in this case is box 0,0. If this box had 
been subdivided, the child container is found by 
shifting the indexes to the left and using the MSBs 
once more. In the particular example shown above, the 
procedure ends after the indexes had been shifted three 
times. The sequence of boxes is then as follows: 00, 
00, 10, 11. 

Most vertices can be eliminated by looking 
only at those inside boxes whose intersections with the 
given shape is not null. To identify these boxes, one 
possibility is to find the smallest box in the OctTree 
index space enclosing the edge plus a neighborhood 
region function of epsilon. To do so, it is first 
necessary to map the edge equation from Cartesian to 
the OctTree index space. The edge is defined as the 
parameter a as: 

(x 0 ,y 0 ,z 0 ) + ot* (Ax, Ay , Az) , with || Ax, Ay, Az|| 2 = 1 

and o £ a < and it starts in P 0 = (x 0 ,y 0 ,z 0 ) (21) 
and ends in P x = (x 0 ,y 0 ,z 0 )+ (Ax,Ay,Az) . 
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Assuming all positive Ax,Ay,Az 4 , the maximum and minimum 
indexes of a box, aligned with the system axis and 
enclosing the edges are: 



P -x 

* 0,jy min 
. ^max ^min 



, int 



( B, - v 

* 0,y snun 
V ^max y min 



■2" 



,int 



p uy-y 



w 3^ max ^min 



2" Lint 



P - v 

■* Q,z ^min 2" 
V 3^ max y min J 



min 



2" Lint 



p>,-y, 



min 



<^max ^min /j 



(22) 
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Since the point must fall in an epsilon- 
neighborhood of the segment, the delta is calculated in 
the OctTree index space as follows: 



Ai = int 



f \ 
- 2" +0.5 



V X max X rmn J 



(23) 



and then subtracted and added to the edge-enclosing 
box. The closest elements will be searched among all 
the vertices belonging to any of the sub-boxes 

10 satisfying the requisite equation: 

^ n -(A/,A/,A*,)<7<;L* + (A/,A*,A/) (24) 
In the example illustrated in Figure 13, this 
would require checking all boxes whose indexes satisfy 
the following equation. 

15 0011,0001-1,1 </ < 0111,0101 + 1,1 => 0010, 0000<T < 1000,0110 < 25 ) 
The sequence of boxes to be checked is 
generated. It is sufficient to traverse, the tree, 
depth-first. For the particular example described 
above, the tree is shown in Figure 14. The search 

2 0 starts browsing the leftmost branch of the 'tree until 
all leaves are traversed, and then continues exploring 
the second leftmost branch of the tree. The last two 
branches are discarded because their indexes are 
outside the allowed range, that is, the following 

25 conditions are not satisfied: 

00 1 0, 0000 < \xxx, Oxxxx < 1 000, 0110 and (26) 

00 1 0, 0000 < \xxx, Ixxxx < 1 000, 0110 (27) 
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A more efficient alternative is to only check 
the boxes that are in the close neighborhood of the 
edge, as shown in Figure 15. The equation of the edge 
in Cartesian space, which is 



P(a)= (x 0 ,y 0 ,z 0 ) + a- (Ax,Ay,Az) 
with || Ax, Ay, Az|| 2 = 1 and 0^a<a w 
becomes the OctTree index space 



(28) 



x - x y n v 



X — X V — V- 2 — Z . 
max rain * max nun max nun 



i\d) = 

r a.. *. > 



Ax Ay Az 



i Z — z. z — z. z — z . j 

\ max mm max nun max mm / 



(29) 



Using this equation, it is now possible to 
traverse the tree and for each box to determine if it 
contains a portion of the edge. Even if the technique 
requires additional computations to determine the list 
of boxes to be included in the search, it culls some of 
15 the boxes that would otherwise have been selected and 
as a consequence reduces the number of vertices whose 
distance from the target needs to be evaluated. On 
average, the overall number of calculations is less 
than in the previous implementation. Since OctTree are 
2 0 three-dimensional binary trees, the searching algorithm 
in accordance with the present invention has a 
logarithmic, base 2, complexity. 

While the invention has been described in 
connection with one or more embodiments, it is to be 
2 5 understood that the specific mechanisms and techniques 
which have been described are merely illustrative of 
the principles of the invention. Numerous 
modifications may be made. to the methods and apparatus 
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described without departing from the spirit and scope 
of the invention as defined by the appended claims. 
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